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0^ ' We apply Chiral Perturbation Theory at one loop to analyze the general formula 

r~| ' for the 'k~^tt~ atom lifetime derived recently in the framework of QCD [1]. The 

C^, corresponding analytic expression is investigated numerically, and compared with 

Q^' recent work in the literature. 
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1. The DIRAC collaboration at CERN [2] aims to measure the lifetime of the 
TT+TT" atom (pionium) in its ground state at the 10% level. This atom decays 
predominantly into two neutral pions, V = r2^o + T2-y + . . . , with T2-y/T2TrO ~ 
4 • 10~^ [2]. The measurement of r2,r0 allows one [1,3-5] to determine the 
difference Oq — 02 of the strong S'-wave vrvr scattering lengths with isospin 
/ = 0, 2. One may then confront the predictions for this quantity obtained in 
standard ChPT [6,7] with the lifetime measurement, and furthermore analyze 
the nature of spontaneous chiral symmetry breaking in QCD [8]. In order to 
perform these investigations, one needs to know the theoretical expression for 
the width of pionium with a precision that properly matches the accuracy of 
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the lifetime measurement of DIRAC It is the aim of the present article to 
provide the necessary numerical setting. 

2. The formation of the atom and its subsequent decay into two neutral pions 
is induced by isospin breaking effects in the underlying theory. In the present 
framework (QCD including photons), these are the electromagnetic interac- 
tions and the mass difference of the up and down quarks. In the following, it 
is useful to count a ~ 1/137 and (tti^ — m^)^ as small parameters of order b. 

More than forty years ago, Deser et al. [3] derived the formula for the width 
of the 7r~p atom at leading order in isospin symmetry breaking. Later in 
Refs. [4,5], this result was adapted to the 7r"*"7r~ atomQ. In particular, it was 
shown that - again at leading order in isospin symmetry breaking effects - the 
width V\^Q of pionium is proportional to the square of the difference oq — a2, 



rL°o = I aYiao - a,f ; p^ = (M^, - M^^o - ^M^.c^) V^. 



At leading order in S, the momentum p* becomes J2Mt^+(Mt^+ — M^o) - this is 
the expression used in [4,5]. We prefer to use Eq. (1), because in this manner, 
one disentangles the kinematical corrections - due to the expansion of the 
square root - from true dynamical ones. 

In our recent article [1], we derived a general expression for the pionium life- 
time, that is valid at leading and next-to-leading order in isospin breaking, 



r2.o = '^aYA\l + K). (2) 

The quantities A and K are expanded in powers of 6. In particular, it has 
been shown in [1] thatj^] 



^=-J^ReAV°° + o(5), (3) 

where Rey4^~°° is calculated as follows. One evaluates the relativistic scatter- 
ing amplitude for the process vr+vr^ -^ 7r°7r° at order S near threshold and 
removes the (divergent) Coulomb phase. The real part of this matrix element 
develops singularities that behave like |p|~^ and In2|p|/M^+ near threshold (p 



^ There are a few misprints in Eq. (6) for the pionium decay rate in Ref. [4]. The 
correct result is displayed in Ref. [5]. 

^ We use throughout the Landau symbols 0(x) [o(x)] for quantities that vanish 
Uke X [faster than x] when x tends to zero. Furthermore, it is understood that this 
holds modulo logarithmic terms, i.e. we write also 0{x) for xlnx. 



denotes the center of mass momentum of the charged pions). The remainder, 
evaluated at the tt+tt" threshold p = 0, equals Rey4j^7°°- It contains terms of 
order 6^ and S and is normalized such that, in the isospin symmetry limit, 
A = aQ — a2. Finally, the quantity K starts at order alna - its explicit ex- 
pression up to and including terms of order 6 is given by 



K = ^^ (ao + 2a2f - y (In a - 1) (2ao + a^) + o{5) , 

A^ = M^-M^o. (4) 



In the derivation of Eqs. (2) - (4), the chiral expansion has not been used - 
details will be provided in a forthcoming publication [9] . 

The following remarks are in order. First, it turns out that numerically, the 
correction factor (1 + K) is very close to unity, see below. The present un- 
certainties in the values of the scattering lengths that enter the expression 
for K are thus of no significance in this respect. Second, it follows that the 
measurement of the pionium lifetime amounts to a precise measurement of the 
quantity A, that is not affected by use of any chiral expansion. 

3. In order to extract tttt scattering lengths from experimental information on 
the width, we invoke ChPT and relate the amplitude A and the difference of 
the S'-wave scattering lengths Cq — 02 order by order in the chiral expansion. 
This may be achieved as follows. First, we expand A in powers of the isospin 
breaking parameter 5. Because the scattering amplitude tt^tt" — ^ ti'^tt^ in 
QCD does not contain terms linear in the quark mass difference m^ — md-, we 
write 



A = a^ ~ a2 + hi {iTbd - niuY + /?.2 « + o{5) , (5) 

where ao and 02 denote, hy definition^ the strong tttt scattering lengths eval- 
uated in QCD at e = 0, ruu = fnd- The quark masses are tuned such that 
the pion mass in the isospin-symmetric world coincides with the charged pion 
mass. The values of other parameters are not changed when performing this 
isospin symmetry limit. For illustration, the difference Oq — 02 at one loop is 
given by [7] 



q i\/t2 / 114-2 _ _ \ 

«o - «2 = ^ (^1 + ^^^ {33 + 8/, + I6/2 - 3/3}] + 0{M!,) , (6) 

where k denotes the running coupling constant /[ at scale /i = M^+, and F 
stands for the pion decay constant in the chiral limit, F ~ 88 MeV [6]. 



Second, the coefficients hi are evaluated in the framework of chiral perturba- 
tion theory. This calculation is straightforward - one only needs to determine 
in the Ti^n' -^ tt^tt'^ amplitude the terms proportional to a and (m^ — niuY- 
Fortunately, the result at order e^p^ is already available in the literature [10]. 
At this order, hi vanishes. Indeed, hi is at least of order rh in any order of the 
chiral expansion, with rh = {rriu + md)/2. The expression for /i2 can be read 
off from amplitude given in Ref. [10]. The result is 



hi = 0{m) , 

h2 = hA + h^ + 0{m^) , 



h^ 
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p{ki) . 



(7) 



Here, p{ki) denotes a particular combination of electromagnetic low-energy 
constants k^ that occur in the effective Lagrangian of SU{2) x SU{2) at order 
e^p^ [10,11]. Using the Lagrangian of Ref. [10], one has 



p{ki) = -30 + 9ki + Qk^ + 2k(i + ks + -Z{ki + 2k2 + Qk^ + 12A;6 - 6A;8) 



and 



A! 



A 



■rT\mu=m^ ; 



87raF2 



mu=m(;=0 



(9) 



The quantities ki denote again the running couplings A;[ at scale /x = M^+. 
Note that according to our counting, the quantity Lt, = \n{M^+/M^o) intro- 
duced in Ref. [10], is of order 6 and hence does not contribute to /i2. 



4. In the numerical analysis, it is convenient to relate the particular combina- 
tion of electromagnetic low-energy constants ki that appears in Eq. (8), to the 
low-energy constants K^ in the SU{3) x SU{3) version of ChPT. Estimates 
for the numerical values of K^ are available in the literature [12-14], whereas 
the ki have not yet been determined to the best of our knowledge. The rela- 
tion between ki and K^ is straightforward to work out: one e.g. evaluates the 
amplitude 7r~^TT~ -^ tt^tt^ and M^o in the framework of SU{3) x SU{3) and ex- 
pands the result at small momenta and small quark masses, p,mu, 171,1 *^ it^s- 
In this limit, the expression for the SU{3) x SU{3) amplitude goes over into 
the one given in [10], provided that one sets 



p{K,) = ^ (-6(ir[ + k;) + mi - 5k; + k^ + q{ki + ki, + ki,] 

- 18 + 28Z ln^i-2Z ln^-^ + 1 -30. 10 

Here, Bq is related to the chiral condensate in SUi^) x ^[/(S) and m^ denotes 
the strange quark mass [6]. Furthermore, we have used that, at this order in the 
low-energy expansion, Z equals its SUi^) x SUi^ analog Zq. The couplings 
Kl can be expressed [13] as a convolution of a QCD correlation function 
with the photon propagator, plus a contribution from QED counterterms. We 
have checked that Pik^ is independent of the QCD scale /xq that must be 
introduced in the QCD Lagrangian after taking into account electromagnetic 
effects [13,14]0. 

Using the relation between F and F,r in QCD [6], 



F. = F(l + ^^r, + 0(m^)), (11) 

the decay width of pionium can finally be expressed in the form 



r2.o = \ c?f (ao - a2 + ef (1 + i^) , (12) 



where 



e = a/iA + Oih. 



"1 



3^e.m. . ^2 

1 + 



327rF2 I 127r2F2 



23 J 3^ 3^ 

h /l + - /s + - ^4 

8 4 2 



The ellipses denote terms of 0(m(m^—m^i)^, am^) and of o((5). We expect these 
terms to give a negligible contribution to e. Note that the term proportional 
to ^4 cancels in the sum a/iA + oih^ at O(e^p^). The formulae (12) and (13) 
will be used below for the numerical analysis. 

5. For the numerical analysis of Eqs. (12), (13), we use the following values 
for the various quantities that occur in these expressions. First, we recall 
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that the non-electromagnetic part of the pion mass difference is tiny, of order 
~ 0.1 MeV [15]. Therefore, we identify A^'™' with the experimentally measured 
total shift A^. Similarly, the value of Z is determined from Eq. (9) using 
the observed value of A^. Further, in the calculations we replace m^-Bo by 
M|-+ — M^+/2, according to our definition of the isospin symmetry limit. 
The values used for the low-energy constants are: F^ = 93.2 MeV [6,16], 
h = -2.3 ± 3.7, k = 6.0 ± 1.3, h = 2.9 ± 2.4 [6], h = 4:A ± 0.3 [17]. For 
K^{fi), we use the values given by Baur and Urech in Ref. [12, Table 1]: K^ = 
-6.4, Kl = 6.4, Kl = -6.2, K^ = 19.9, K^ = 8.6, K^ = K{q = 0, K[^ = 0.6 
(in units of 10"^). We evaluate P{K.j) at scale /x = Mp. Further, we attribute 
an uncertainty 2/167r^ - that stems from dimensional arguments - to each 
Kl. The values of Kl obtained both by Moussallam [13] and by Bijnens and 
Prades [14] , lie then within the uncertainties attributed. In the final expression 
for e, the uncertainties coming from /j and Kl, are added quadratically. Finally, 
we use for a^ and a2 the values corresponding to set 2 in Ref. [16]: ag = 0.206, 
a2 = —0.0443, without attributing any error. These scattering lengths enter 
only the correction K in (12), which is very small. 

Inserting all this in e, and omitting the terms of higher order in the chiral 
expansion and in isospin violating effects [indicated by the ellipses in (13)], we 
arrive at 



(0.58 ±0.16) -lo- 



ir = 1.07-10" 



(14) 



The error in this result includes only the uncertainty in the values of the 
low-energy constants Kl and /,. Resonance saturation introduces a scale de- 
pendence in the final result for e. If saturation is assumed at ;U = 500 MeV 
(/i = 1 GeV), we find e = 0.51 - lO'^ (e = 0.62 - lO'^). Eqs. (12) - (14) are 
the main result of this letter. We add the following remarks concerning this 
analysis. 

i) The term ah/\ {ah^) contributes with 0.51 - 10~^ (0.07 - 10~^) to e, see 
Eq. (13). Therefore, ah/^ is by far the dominant effect in the improvement of 
the leading order formula (1). The dominant contribution in K stems from 
the logarithmic term ~ a In a. 

ii) At one-loop order, ah^^ can be written in the following form. 



ah A 



A! 



3M2 



ao — a2 + 



Mt 



2567r3Fi 



9 + 4(^1 - h) + y Is 



(15) 



For the values of the low-energy constants /i, I2 and ^3 given above, the second 
term in the curly bracket turns out to be very small, of order —6 - 10~^. 
Consequently, the one-loop correction to the ratio ah^/i^ao — 02) is tiny, and 
the latter is dominated by the tree-level contribution A^/3M^+. 



iii) If one is willing to rely on the numerical values of oq and a2 displayed above, 
one may evaluate the lifetime of pionium in the ground state: r = 3.25-10"^^ s. 
The correction to the lowest-order formula (1) by Deser et al. then becomes 
(r2^o — r2^))/r2^o = O.OSS. We prefer not to attach any error to these numbers, 
because they are based on the numerical values of oq and 02, and these are 
not yet known with sufficient accuracy. 

6. We now compare our result with recent work in the literature and start 
the discussion with Refs. [18], where hadronic atoms have been studied in the 
framework of potential scattering theory. The correction to the leading-order 
formula for the width has been worked out in [18] numerically - a comparison 
with our analytic result is thus not possible. The numerical result quoted in [18] 
differs significantly from ours: isospin violating corrections to the leading order 
result increase the lifetime according to these authors, in contrast to Eqs. (12), 
(14). The discrepancy does not come as a surprise - in the present form, the 
potential approach [18] does not reproduce all isospin-breaking terms that are 
present in the Standard Model. The leading correction in ChPT stems from 
tuning the quark mass in the isospin-symmetric phase such that the pion mass 
in the isospin-symmetric world coincides with the charged pion mass. Since 
the pion-pion interaction depends on the quark mass, this effectively leads to 
a change in the tttt potential when the isospin limit is considered. In addition, 
the potential model in [18] does not take into account the direct quark-photon 
effects encoded in the low-energy constants ki in our approach. We believe 
that, for a consistent calculation of the corrections to the tt'^tt^ atom decay 
width, the potential in the scattering theory approach should be matched to 
ChPT in the isospin violating phase. This procedure would guarantee that the 
above mentioned effects are included. 

7. In order to present a coherent comparison with other calculations performed 
in the framework of Quantum Field Theory, it is useful to expand also the 
quantity K, similarly to A in (5), 

K = fi {nid - muY + f20i\na + fsa + o(S), (16) 

and to rewrite the decay width (2) in the form 

R = ao - a2 + Via In a + r2a + rs{md - niuf + o{5) , (17) 

where the coefficients in the expansion are given by 



ri = -{aQ-a2)f2^ 



r2 = h2 + -(ao-a2)/3, 

r3 = hi + -{ao-a2)fi. (18) 

The comparison with other approaches becomes then rather easy: one com- 
pares the analytic expressions for Tj at a given order in the chiral expansion. 

The correction to the pionium decay width were evaluated in Refs. [19] by 
use of 3D constraint theory equations, and in Refs. [20] in the Bethe-Salpeter 
approach. The term proportional to a^/A^r that emerges from the expansion 
of p*, was omitted in the final expressions for the decay width in these papers. 
Note that though this term is algebraically of order 6, the numerical effect 
coming from it is negligible. Further, rewriting the expressions found in these 
works in the form (17,18), we find that the coefficient ri coincides to all orders 
in the chiral expansion with our result, whereas r2 coincides with the above 
result up to and including terms of order p^. The matching relation (10) does 
however not agree with the one proposed in Refs. [19]. The expressions in 
Refs. [19] and [20] contain some of the higher-order terms in the quark mass 
expansion and in isospin breaking effects. These cannot be reliably predicted 
based on O(e^p^) calculations alone. The effect of those terms is however 
small. The numerical values of the corrections in [19,20] almost coincide with 
our result. 

Labelle and Buckley [21] were the first to apply the non relativistic effective 
Lagrangian approach to the problem of hadronic atoms. They evaluated the 
leading order term, and in addition the correction due to electron vacuum 
polarization. This contribution is analytically of order a^ and thus beyond the 
accuracy considered here. Although this correction is potentially large, because 
it is proportional to {Mt^+ /meiectron)<^'^ — 2q;, it amounts to a tiny contribution 
3.9-10""^ to e, which one may safely ignore. The vacuum polarization correction 
to the 7r^7r~ atom decay width was obtained independently in Ref. [22] - 
the result agrees to that of [21]. Note that the part of the correction due to 
vacuum polarization, with no Coulomb corrections in the intermediate state, 
was considered also in Refs. [20] , and the result agrees with the corresponding 
one given in Ref. [21]. 

Recently, Kong and Ravndal [23,24] and Holstein [26] have applied non rela- 
tivistic effective Lagrangian techniques to this problem as well. In Refs. [23,24], 
Coulomb corrections are not taken into account, which amounts to ri = 0. The 
relativistic correction found in Ref. [24] amounts to rewriting the phase space 



factor w2Mjr+(M^+ — Mj^o) as JM"^+ — M^o - this term is thus included in our 
approach. On the other hand, we do not agree in the value of the coefficient 
r2 given in Ref. [23] already at leading order in the quark mass expansion. 
The discrepancy is due to the fact that - in [23] - the matching of the ef- 



fective couplings in the non relativistic Lagrangian has not been carried out 
with a precision that is required to pin down all terms at this order in the 
expression for the width. If a complete matching at 0{S) is performed, agree- 
ment is achieved at the leading order in quark mass expansion [25] . The result 
of Ref. [26] for r2 agrees with Ref. [23] - hence, it also shares its shortcom- 
ings. The contribution corresponding to ri is omitted in the final expression 
for the width in [26, Eq.(95)], although at previous stages in that work, the 
Coulomb corrections had been discussed. Finally, effects from rriu 7^ rriii are 
not disentangled in [23,24,26]. 

In Ref. [27], an effective non relativistic Lagrangian was used to derive the ex- 
pression for the decay width in terms of effective couplings in that Lagrangian. 
The result agrees with ours [1, Eq.(ll)]. The matching to the relativistic am- 
plitude is carried out in [27] order by order in the chiral expansion. We agree 
with the result of [27] for the coefficient ri at leading order in the chiral ex- 
pansion. The result for the coefficient r2 is not explicitly given in [27], and a 
comparison is therefore not possible. 

The result presented in Ref. [28] differs from ours in many respect. In particu- 
lar, the expression for the decay width given in Ref. [28] contains an ultraviolet 
divergence that is regularized by introducing an explicit cutoff. A systematic 
renormalization procedure is not discussed. For this reason, the result of [28] 
can not be compared with the present one. 

Finally, we comment on the work of Ref. [29]. First we note that in our ap- 
proach, the decay width is calculated by considering the elastic 7r°7r° — > 7r°7r° 
amplitude as a function of the energy E of the ii^ii^ pair. This amplitude 
develops a pole at E = Ep on the second Riemann sheet - the imaginary part 
Iin{Ep) is related to the width in the standard manner. As far as we can see, 
the approaches proposed in [23,24,26] amount to the same definition of the 
width, because it is identified in these works with the imaginary part of the 
level shift in Rayleigh-Schrodinger perturbation theory. On the other hand, 
Ref. [29] amounts to an alternative definition of r27r0 - a direct comparison 
with the present work is therefore not possible. 

8. In conclusion, we have analyzed analytically and numerically the general 
formula [1] for the decay width of the pionium ground state at one-loop order 
in ChPT, and we have compared our result with other work available in the 
literature. We have in particular identified the reason for discrepancies of our 
result with other approaches. As we have shown, a precise determination of the 
pionium lifetime allows one to measure the amplitude A in Eq. (2), where no 
chiral expansion is used. By use of ChPT, one may determine the combination 
do — 02 of S'-wave scattering lengths from Eqs. (12) - (14), that constitute the 
main result of this letter. As we expect the corrections due to higher orders 
in the chiral expansion to be small, we infer from this result that an accurate 
determination of Cq — a2 from a precise lifetime measurement is indeed feasible. 
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